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A method of and an apparatus for processing seismic data 
Field of the Invention 

The present invention relates to a method of processing seismic data, in particular 
nmlticomponent seismic data, to remove unwanted events from the acquired data. It 
also relates to an apparatus for processing seismic data. It also relates to a method of 
and an apparatus for calculating a de-signature and de-multiple operator or a de- 
signature operator. 

The Related Art 

Figure 1(a) is a schematic illustration of the principles of a seismic survey. This seismic 

* 

survey is intended to provide information about a target geological reflector 3 disposed 
within the earth's interior. 

The seismic survey shown in Figure 1 includes a seismic source 4 disposed on the 
earth's surface 1. A seismic sensor 5, referred to hereinafter as a "receiver 5 ', is also 
disposed on the earth's surface, at a distance from the seismic source 4. In use, the 
seismic source 4 is caused to actuate a pulse of seismic energy, and emitted seismic 
energy is detected by the receiver 5. 

Figure 1(a) illustrates a land-based seismic survey. However, seismic surveying of this 
general type is not restricted to land, and may be carried out in a marine environment or 
in the land-sea transition zone. For example, marine seismic surveying arrangements 
are known in which one or more seismic sources are towed by a survey vessel; in such 
an arrangement, the receivers maybe disposed on the sea-bed (a so-called "ocean 
bottom cable" survey or OBC survey), or the receivers may also be towed by a survey 
vessel. Furthermore, a land-based seismic survey is not limited to the configuration 
shown in Figure 1(a), and it is possible for the seismic source or the seismic receiver to 
be disposed within the earth's surface. For example, in a vertical seismic profile (VSP) 
seismic survey a seismic source is placed on the earth's surface and a receiver is 
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disposed within the earth's interior in a bore hole. In a reverse VSP survey a seismic 
source is disposed within a bore hole, and a receiver is disposed on the earth's surface. 

Only one seismic source 4 and one seismic receiver 5 are shown in Figure 1(a) for ease 
of explanation but, in general, a practical seismic survey will contain an array of sources 
and an array of receivers. 

One problem that arises in seismic surveying is that seismic energy may travel from the 
source to the receiver along many paths. One reason for this is that many other 
reflectors exist within the earth in addition to the target reflector. In Figure 1(a) this is 
illustrated schematically by a reflector 2 that overlies the target reflector 3. These 
additional reflectors generate paths of seismic energy from the source to the receiver 
that involve a reflection at a reflector other than the target reflector. Another reason for 
the existence of multiple paths of seismic energy is that seismic energy that is 
propagating upwardly within the earth will undergo reflection at the earth's surface 1, 
owing to the different seismic properties of the earth and the air. This leads to the 
existence of paths of seismic energy from the source to the receiver that involve more 
than a single reflection at the target reflector. These paths give rise to unwanted events 
in seismic data acquired at the receiver. An event in acquired seismic data that relates 
to seismic energy that has undergone multiple reflections will hereinafter be referred to 
a "multiple event". 

Figure 1(a) illustrates the primary path of seismic energy for this surveying 
arrangement, in which the path of seismic energy from the source 4 to the receiver 5 
involves only a single reflection, at the target reflector 3. (Refraction at the overlying 
reflector 2 has been omitted from Figure 1 (a) for clarity.) Although seismic energy 
propagating along the primary path passes through the overlying reflector 2 on its 
downwards path from the source 4 to the target reflector 3, and again on its upwards 
path from the target reflector 3 to the receiver 5, the primary path does not involve 
reflection at the overlying reflector 2. In an ideal seismic survey, only seismic energy 
that travelled along the primary path would be detected by the receiver 5. 
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In a practical seismic survey, the receiver 5 will detect seismic energy that has travelled 
from the source 4 along many paths other than the primary path. Examples of these 
other paths of seismic energy are shown in Figure 1(b) to 1(d). Figures 1(b) to 1(d) 
illustrate paths of seismic energy that involve more than one reflection, and these are 
known as "multiple events". In Figure 1(b) downwardly propagating seismic energy 
from the source 4 is reflected by the overlying reflector 2 so that it travels upwards to 
the earth's surface 1. The seismic energy is further reflected downwards at the earth's 
surface, and is then incident on the target reflector. An event of this general type is 
known as a "source-leg multiple", since the additional reflections occur in the path of 
seismic energy from the source to the target reflector. 

Figure 1(c) illustrates a seismic energy path in which the seismic energy travels direct 
from the source to the target reflector 3, and is reflected upwards at the target reflector 
3. However, the reflected seismic energy is not incident direct on the receiver, but is 
reflected downwards at the earth's surface, and is then reflected upwards at the 
overlying reflector 2 before reaching the receiver. A seismic path of this type is known 
as a "receiver-leg multiple", since the additional reflections occur on the path of seismic 
energy from the target reflector to the receiver. 

Figure 1(d) illustrates a seismic path in which seismic energy from the source is 
incident on the target reflector 3, is reflected upwards to the earth's surface, is reflected 
downwardly and undergoes a further reflection at the target reflector 3 before reaching 
the receiver. In this seismic energy path, the additional paths occur between the path of 
seismic energy from the source to the target reflector and the path of energy from the 
target reflector to the receiver. 

The seismic energy acquired at the receiver in a practical seismic survey will include 
events corresponding to the desired primary path 1(a), but will also contain events 
relating to unwanted multiple paths such as the paths shown in Figures 1(b) to 1(d). In 
order to provide accurate information about the target reflector it is desirable to be able 
to identify and remove multiple events from the seismic energy acquired at the receiver. 
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The seismic sources used in a land-based survey are normally vibrated or explosive 
sources. If vibrators are used it is possible to perform a multi-component survey, using 
a multi-component vibrator that produces three orthogonal source motions (two in 
orthogonal horizontal directions and one in the vertical direction). If a seismic receiver 
is used that can record particle motion in three orthogonal directions, then it is possible 
to perform a 3C x 3C (or 9C) seismic survey. A suitable receiver for this is one that can 
measure three orthogonal components of the particle motion at the receiver, for example 
a receiver that contains three orthogonal geophones - two geophones for measuring two 
orthogonal horizontal components of the particle motion at the receiver, and a third 
geophone for measuring the vertical component of the particle motion at the receiver, 

A further problem in analysing the results of a seismic survey is that multi-component 
vibrators, which operate by imposing tractions on the earth's surface, emit three distinct 
wave types, known as P-, Sv- and Sh-waves (P-waves are pressure waves, and Sv and 
Sh waves are shear waves). The relative amplitudes of these different wave-types in 
seismic energy emitted by a multi-component vibrator vary depending on the direction 
of propagation of the seismic energy. A multi-component receiver records the three 
wave types, with a sensitivity that depends on the angle of incidence of the received 
seismic energy. When a geophone measures a component of the wavefield at the 
earth's surface, both P-waves and the two-types of S-waves are recorded without 
distinction. 

This is illustrated schematically in Figure 4(a). Figure 4(a) shows a land-based seismic 
survey in which a three-component seismic source 4 - in this case a multi-component 
vibrator - and a three-component seismic receiver 5 are disposed on the earth's surface. 
As indicated in Figure 4(a), the seismic source 4 emits both P waves and S waves (only 
one S-wave type is shown for clarity), and the receiver 5 detects both P-waves and S- 
waves. The four views in Figure 4(a) show, from left to right: 

(i) Seismic energy generated by a horizontal source motion at the vibrator 4, and 
being received by a horizontally-oriented geophone at the receiver 5; 
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(ii) Seismic energy being generated by a horizontal source motion at the vibrator 4, 
and being detected by a vertically— oriented geophone at the receiver 5; 

(iii) Seismic energy being generated by a vertical source motion at the vibrator 4 and 
being detected by a horizontally— oriented geophone at the receiver 5; and 

(iv) Seismic energy being generated by a vertical source motion at the vibrator 4 and 
being detected by a vertically-oriented geophone at the receiver 5. 

Figure 4(a) illustrates only one horizontal component and one vertical component of the 
source motion at the vibrator 4 and one horizontal geophone component and one 
vertical geophone component at the receiver 5. As noted above, in a full multi- 
component survey the vibrator 4 will further generate seismic energy by a source 
motion out of the plane of the paper, and the receiver will further comprise a third 
geophone that detects the particle motion along a line out of the plane of the paper. In 
total, thus, there are 9 combinations produced by the three orthogonal source motions 
generated at the source and the three orthogonal geophones at the receiver. 

As noted above, many seismic receivers record P- and S-waves without distinction, so 
that a seismic trace acquired at the receiver will includes events due to received P-waves 
and events due to received S-waves. In many cases it is desirable to separate the P- 
events in a seismic trace from the S -events, since this provides additional information 
about the earth's interior. In many cases, a geological structure will have a different 
effect on P-waves than on S-waves. The process of separating the P-events in a seismic 
trace from the S-events is generally referred to as decomposing the seismic trace into its 
P- and S-components. 

There are a number of prior art approaches to decomposing seismic energy acquired at a 
receiver into P-components and S-components. There are also a number of prior 
approaches to eliminating the effect of multiple reflections from the acquired seismic 
data. 
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For land-seismic data, CP. A. Wapenaar et al, in "Decomposition of multicomponent 
seismic data into primary P- and S-wave responses", Geophys. Prosp. Vol. 38, pp633- 
661 (1990) and P. Herrmann, in "Decomposition of multicomponent measurements into 
P- and S-waves", Ph.D. thesis, Delft University of Technology (1992) have given an 
elastic decomposition scheme which decomposes data both at the receiver side 
(common shot gather) and at the source side (common receiver gather). The wave field 
decomposition applied to a common source gather along the receiver line replaces the 
original particle velocity detectors by pure P- and S- wave detectors. The wave field 
decomposition applied to a common receiver gather similarly replaces the original 
vibrator sources by pure P- and S- wave sources. Thus, the total decomposition scheme 
provides 9 different data sets: P, Sv and Sh processed gathers from simulated P, Sv and 
Sh sources. 

After the wave field decomposition, Wapenaar et al (1990, supra) apply an inverse 
scheme that eliminates the response of the earth's surface from the decomposed results. 
In this technique, the earth's surface is assumed to be an interface between a solid and a 
vacuum. Eliminating the response of the earth's surface eliminates all multiple paths 
that involve reflection at the earth's surface from the decomposed data, so giving the 
primary P-and S- wave responses of the earth's interior. 

This prior technique requires that the seismic sources are point sources and have a 
known wavelet or have a wavelet that can be estimated from the acquired data. 

The above prior art technique has been generalised by E. Holvik and L. Amundsen, in 
"decomposition of multi component sea floor data into primary PP, PS, SP, and SS 
wave responses", Expanded Abstracts of 68 th Annual frit Mtg of Society of Exploration 
Geophysicists, pp2040 - 2043 (1998), to the case of a marine seismic survey in which 
ideal vibrators (traction sources) and geophones are deployed at the sea floor. This 
technique relates to decomposing the acquired data, and removing events' arising to 
reflections within the water layer. 
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A further prior art technique has been developed by K. Matson and A. Weglein, in 
* 'Removing of elastic interface multiples from land and ocean bottom seismic data using 
inverse scattering" in Expanded Abstracts of 66th Annual Int Mtg of Society of 
Acceleration of Geophysicists, ppl526-1529 (1996), and by K. Matson in "An inverse 
scattering series method for attenuating elastic multiples from multicomponent land and 
ocean bottom seismic data", Ph.D. thesis, University of British Columbia (1997). In 
this technique, inverse scattering series are used to develop elastic schemes that 
attenuate multiple reflections from multi-component land or ocean bottom seismic data. 

The techniques developed by Holvik and Amundsen and by Matson and Weglein again 
require that the seismic sources must be point sources with a wavelet that is known or 
that can be estimated from the acquired data. 

L. Amundsen has proposed, in "Geophysics" Vol 66 pp 327-341 (2001), a method of 
eliminating free-surface multiples from marine seismic data acquired using a 
multicomponent seismic receiver disposed in the water column or on the sea-bed. 
However, this method does not remove the effects of multiple reflections associated 
with the sea-bed. 

Summary of the Invention 

A first aspect of the present invention provides a method of processing multi- 
component seismic data, the data having been acquired by emitting multi-component 
seismic energy at a source location; and acquiring seismic data at a multi-component 
seismic receiver located at a greater depth than the source location, the method 
comprising the steps of: decomposing the acquired seismic data into up-going 
constituents and down-going constituents; and calculating a de-signature and de- 
multiple operator from the down- going constituents of the acquired seismic data and 
from properties of the medium surrounding the receiver. 

The de-signature and de-multiple operator of the invention is effective at attenuating or 
completely removing the effects of the overburden from the seismic data. It is effective 
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at removing from the seismic data all multiple reflections associated with any interface 
above the level of the receiver or with any interface at the receiver level. (The term 
"interface" covers any discontinuity in acoustic or elastic properties that causes partial 
reflection of seismic energy.) 

It is also effective at attenuating or completely removing the effects of the source 
radiation characteristics (or the source "signature") from the data. 

The method may comprise processing the acquired seismic data using the de-signature 
and de-multiple operator thereby to attenuate or remove seismic events arising from 
multiple reflections. This maybe done by processing the up-going constituents of the 
acquired seismic data using the de-signature and de-multiple operator. Alternatively, 
the de-signature and de-multiple operator may be applied to the complete (that is, 
undecomposed) acquired seismic data. 

The step of processing the acquired seismic data may comprise selecting a desired 
seismic signature for the source. 

The method may further comprise decomposing the seismic data into P-wave and/or S- 
wave data. This decomposition step may be performed on the de-signatured, de- 
multtpled data. Alternatively, the acquired data may be decomposed into P-wave and/or 
S-wave data before the de-signature and de-multiple operator is applied. 

The decomposition into P-wave and/or S-wave data may be a receiver-side 
decomposition and/or a source-side decomposition. 

A second aspect of the invention provides a method of processing multi-component 
seismic data, the data having been acquired by emitting multi-component seismic 
energy at a source location; and acquiring seismic data at a multi-component seismic 
receiver located at a greater depth than the source location, the method comprising the 
steps of: 
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decomposing the acquired seismic data into up-going constituents and down- 
going constituents; and 

calculating a de-signature operator from the initial down-going constituents of 
the acquired seismic data and from properties of the medium surrounding the receiver. 

The initial down-going constituents of the seismic data acquired at the receiver relate to 
seismic energy that has travelled direct from the source to the receiver without 
undergoing any reflection. They may therefore be used to calculate a de-signature 
operator for attenuating or eliminating the effects of the radiation characteristics of the 
seismic source on the acquired data. 

The methods of the invention may be applied to pre-existing seismic data. 
Alternatively, the method may further comprises the steps of: emitting multi-component 
seismic energy at a source location; and acquiring seismic data at a multi-component 
seismic receiver located at a greater depth than the source location. 

Other aspects of the invention are claimed in independent claims 10, 12, 13, 14, 16, 23, 
24, 26, 27,28 and 30. 

Brief Description of the Drawings 

Preferred embodiments of the present invention will now be described by way of 
illustrative example with reference to the accompanying figures in which: 

Figure 1(a) is a schematic view of a seismic survey showing the desired path of seismic 
energy; 

Figures 1(b), 1(c) and 1(d) are schematic views of a seismic survey showing three 
undesired paths of seismic energy; 
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Figure 2 is a flow diagram showing possible steps of a method of obtaining a de- 
signature and de-multiple operator according to an embodiment of the present 
invention; 

Figure 3(a) is a schematic view of a real seismic survey; 

Figure 3(b) shows a hypothetical seismic survey corresponding to the seismic survey of 
Figure 3(a); 

Figure 3(c) shows a further hypothetical seismic survey that is a reciprocal of the 
hypothetical seismic survey shown in Figure 3(b); 

Figures 4(a) to 4(d) illustrate the effect of applying a method of the present invention to 
a land seismic survey; 

Figures 5(a) and 5(b) illustrates the effects of applying the method of the present 
invention to an ocean bottom seismic survey; 

Figures 6(a) and 6(b) illustrate the effect of applying a method of the present invention 
to a land seismic survey in which the receiver is deployed within the earth's interior; 

Figures 7(a) to 7(d) show raw synthetic seismic data; 

Figures 8(a) to 8(d) show the up-going component of the synthetic seismic data of 
Figures 7(a) to 7(d); 

Figures 9(a) to 9(d) show the down- going part of the synthetic seismic data of Figures 
7(a) to 7(d); 



Figure 10(a) to 10(d) show the particle velocities obtained from the synthetic seismic 
data of Figures 7(a) to 7(d) using a method of the present invention; 
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Figure 11(a) to 11(d) show thePP, SP, PS and SS wavellelds obtained from the 
synthetic seismic data of Figure 7(a) to 7(d) using a method of the present invention; 

Figure 12 is block schematic view of an apparatus according to the invention. 

Brief Description of the Preferred Embodiments 

The principle of one method of obtaining a de-signature and de-multiple operator of the 
present invention is illustrated in Figures 3(a) to 3(c). Figure 3(a) illustrates an actual 
seismic survey. This seismic survey comprises a receiver x r (14) that is located at a 
depth Zr. The medium within which the seismic survey is carried out is classified into 
the overburden for all depths less than the receiver depth (that is for all z < z r ), and as 
the sub-surface for depths greater than the depth of the receiver (that is, for all z > Zr). 
The overburden is an arbitrary inhomogenous medium. 

A seismic source x s (13) is disposed in the overburden, at a depth less than the depth of 
the receiver. The seismic source emits both up-going and down-going seismic energy. 

The inhomogenous nature of the overburden is indicated schematically by the presence 
of interfaces 10, 11, 12 within the overburden. The seismic properties of the 
overburden change at these interfaces. In the case of a marine seismic survey, the upper 
interface 12 might be the sea-surface (that is, an air/water interface). Because of the 
existence of these interfaces in the overburden, up-going seismic energy is reflected 
downwards in the overburden so that multiple events occur in the seismic data acquired 
at the receiver 14. 

The seismic source x s (13) emits both S- and P- waves. The receiver x r (14) records 
receive S- and P-waves without distinction. The receiver may be any device that 
measures particle motion, such as particle displacement, particle velocity (the time 
derivative of particle displacement), or particle acceleration (the time derivative of 
particle acceleration). One commonly-used seismic receiver is a geophone, but the 
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invention is not limited to processing seismic data acquired using one or more 
geophones as the receiver(s). 

Figure 3(b) illustrates an ideal seismic surveying arrangement that corresponds to 
Figure 3(a). In the idealised arrangement of Figure 3(b), the overburden is an 
homogenous medium, so that up-going seismic energy from the seismic source 
continues to propagate upwards and is not reflected downwards, so that seismic energy 
acquired at the receiver does not contain multiple events. The interfaces 10, 1 1 and 12 
in the overburden of Figure 3(a) have been illustrated in broken lines in Figure 3(b) to 
emphasise that the overburden is homogenous in Figure 3(b). 

In the ideal seismic surveying arrangement of Figure 3(b) the seismic source 13 is 
centred at a location x r ", which is at a depth z/. The depth of the source, Zr", is at an 
infinitesimal distance E above the depth z r at which the receiver 14 is positioned. The 
distance E is vanishingly small. 

« 

Figure 3(c) shows an alternative idealised seismic surveying arrangement. This 
corresponds to the idealised seismic surveying arrangement of Figure 3(b), except that 
the positions of the seismic source 13 and the seismic receiver 14 have been 
interchanged. That is, in Figure 3(c) the seismic source 13 is located at position x r , at 
depth z r , and the seismic receiver 14 is located at location x r " at a depth Zr" which is a 
distance E above the depth Zp of the source. 

Since the seismic surveying arrangement of Figure 3(c) is identical to the seismic 
surveying arrangement of Figure 3(b) except that the positions of the source and 
receiver have been interchanged, the theorem of reciprocity requires that the seismic 
surveying arrangement shown in Figure 3(b) will produce the same seismic data as the 
seismic surveying arrangement shown in Figure 3(c). 

In one embodiment of the invention, Betti's reciprocity theorem is used to transform 
actual seismic data required in a real seismic survey of the type shown in Figure 3(a) 
into the seismic data that would be obtained in an idealised seismic survey of the type 
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shown in figure 3(b). This enables the effect of the inhomogenous overburden to be 
eliminated, thereby eliminating multiple events. 

The present invention also provides decomposition operators, that decompose the 
acquired seismic data into PP, PS, SP, and SS components. 

Bettys reciprocity theorem gives an integral equation relationship between two 
independent elastic wavefields that are defined in a specific volume enclosed by a 
mathematical or physical surface. The relationship between the two wavefields is 
governed by possible differences in the parameters of the medium, by possible 
differences in the source distribution, and by possible differences in boundary 
conditions. In the invention, Betti's theorem is used to provide an integral equation 
procedure that can be used to transform a wavefield that was recorded in an actual 
seismic survey with an overburden response (that is, a survey in which multiple 
reflections occur owing to an inhomogenous overburden) into wavefields that would 
have been recorded in an idealised hypothetical seismic experiment with the overburden 
response absence (the overburden in the hypothetical seismic surveying arrangement is 
homogenous and so produces no response). The present invention makes it possible to 
remove undesired events in the seismic data that arise owing to multiple reflections of 
the types shown in Figures 1(b) to 1(d), so that a multiple-free image can be obtained. 
The present invention enables all multiple reflections associated with any interface 
above or at the level of the receiver to be removed from the seismic data. 

Other than the position of the source, or sources, the present invention does not require 
any source characteristics for the calculation to eliminate the waves scattered from the 
overburden. The radiation characteristics of the physical seismic source(s) is eliminated 
in the transformation from the physical seismic surveying arrangement to the idealised 
hypothetical seismic surveying arrangement. The source in the hypothetical idealised 
surveying arrangement is assumed to be a single point source that is generates 
components orthogonally oriented with respect to one another, and has a desired 
wavelet. 
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This embodiment of the method of eliminating the response of the overburden and the 
effect of the multi-component source radiation pattern from the acquired seismic data is 
referred to hereinafter as "Betti designature and elastic demultiple". A Betti designature 
and elastic demultiple method of the invention has the following advantages: 

(a) It preserves the amplitude of the primary reflection while eliminating all events in 
the seismic data that arise owing to scattering of up-going energy in the 
overburden; 

(b) It requires no knowledge of the medium below the depth of the receiver; 

(c) It requires no knowledge of the medium above the depth of the receiver; and 

(d) It requires information only about the local density and local velocities of elastic 
wave propagation at the depth of the receiver. 

In a method of the present invention recorded seismic data are transformed to new, 
idealised data that would be recorded in a hypothetical seismic survey in which the 
overburden response was absent (that is in which the overburden was completely 
homogenous). The source in this hypothetical seismic survey is a point source, with 
some desired signature radiation characteristics (wavelet). In the processed data, the 
effect of the physical source and its radiation characteristics is absent. In other words, 
the new data have been designatured. 

A method of the present invention may be applied to data obtained in any seismic 
survey that uses a multi-component seismic source and a multi-component seismic 
receiver. Three particular seismic surveying arrangements to which the invention can 
be applied are as follows: 

A. A 3C x 3C (or 9C) land seismic survey, in which the seismic source generates 
three orthogonal source motions (usually two horizontal source motions and one vertical 
source motion). The seismic wavefield is recorded by a receiver that contains three 
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orthogonal geophones, deployed at or below the earth's surface. The three orthogonal 
geophones measure three orthogonal components of the particle velocity vector (usually 
the x-, y- and z-components). 

B. A 3C x 4C (or 12C) marine seismic survey, in which the seismic source generates 
three orthogonal source motions (usually two horizontal source motions and one vertical 
source motion). The seismic wavefield is recorded by a receiver that contains three 
orthogonal geophones, deployed just below the seabed and one hydrophone deployed 
just above the seabed. The three orthogonal geophones measure three orthogonal 
components of the particle velocity vector, and the hydrophone measures the pressure 
field (which is a scalar quantity). 

C. A 3C x 6C (18C) bore hole seismic experiment in which the seismic source again 
generates three separate orthogonal source motions. The source may be disposed on 
land, on the sea floor, or in the water column. The seismic receiver preferably measures 
three orthogonal components of the particle velocity, and three orthogonal comppnents 
of the vertical traction vector or estimates thereof. 

According to the invention, the de-multiple and de-signature operator is obtained from 
the down-going constituents of the acquired seismic data. This means that the 
measurements made at the receiver must be sufficient to allow the acquired wavefield to 
be decomposed at the receiver side into up-going and down-going wave constituents at 
the receiver. In the most general case, a proper wave field decomposition requires that 
the receiver measures the three components of the particle velocity vector as well as the 
three components of the vertical traction vector. For land seismic data, however, the 
vertical traction vector is zero, so that only the three components of the particle velocity 
vector need to be measured. 

For ocean bottom seismic data, the two horizontal components of the vertical traction 
vector are zero, and the vertical component of the vertical traction vector just below the 
sea-bed is equal in magnitude and opposite in sign to the pressure field just above the 
sea floor. Thus, providing a hydrophone just above the sea-bed enables the vertical 
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component of the vertical traction vector (which is the only non-zero component of the 
traction vector) to be determined. 

Applying the present invention to a bore hole seismic survey requires, in principle, that 
the three components of the vertical traction vector must be measured or otherwise be 
known. In most cases, however, the three components of the vertical traction vector are 
not known. In order to overcome this, it is usual to record the three-component particle 
velocity at two different depths within the bore hole. Thus, the receiver system is 
essentially a six-component receiver system, so that a bore hole seismic survey may be 
considered to be an 1 8C survey. 

Figure 2 is a flow diagram showing the principal steps of a method according to the 
present invention. 

The invention may be carried out on suitable pre-existing seismic data. In this case, the 
method would commence with the step (not shown) of retrieving suitable seismic data 
from storage. Alternatively, the method may commence with the step (not shown) of 
performing a multi-component source, multi-component receiver seismic survey to 
acquire suitable seismic data. 

Initially at step SI the local elastic properties of the medium in the vicinity of the 
receiver are determined. For a marine seismic survey, as one example, this step requires 
a determination of the local elastic properties of the sea-bed in the vicinity of the 
receiver. These properties may be determined from measurements, or they may be 
estimated from knowledge of the geological structure of the survey site. 

At step S2, the multi-component seismic data that has been acquired or retrieved from 
storage is decomposed into up-going and down-going wavefield constituents. 
In principle step S2 could be carried out before, or simultaneously with, step SI. It is 
generally preferable to perform step SI before step S2, however, since the local elastic 
properties obtained from step SI may then be used to decompose the seismic data into 
up-going and down- going waves. 



WO 03/036331 



PCT/IB02/05056 



M&C Folio No P51685WO 17 

At step S3 the de-signature and de-multiple operator is calculated from the multi- 
component data and from the elastic properties of the medium. 

At step S4 a desired signature for the seismic source is selected. This desired source 
signature, the multi -component data obtained at S2, and the operator calculated at S3 
are then used to solve an integral equation to find the de-signatured and de-multipled 
wavefield, at step S5. In essence, S5 consists of transforming the seismic data obtained 
in a real seismic survey of the type shown in Figure 3(a) into the data that would be 
obtained in an idealised seismic survey shown in Figure 3(b). In a preferred 
embodiment Betti's theorem is used to transform the actual, measured wavefield 
enclosed by a surface S of the seismic surveying arrangement of Figure 3(a) into the 
wavefield enclosed by the same surface S in the idealised seismic surveying 
arrangement of Figure 3(b). The surface S is made up of the surfaces 2 and Sr shown 
in Figures 3(a) to 3(c). The surface 2 is a horizontally plane surface that is located at 
depth Zr" - that is, located an infinitesimal distance above the depth of the receiver 
(which is located at a depth Zr). The surface Sr is a hemispherical surface having a 
radius of R. 

In Step S5 the de-signature and de-multiple operator may be applied to the up-going 
wavefield constituents obtained at step S 1 . Alternatively, the de-signature and de- 
multiple operator maybe applied to the seismic data as originally- acquired at the 
receiver. 

Finally, at step S6, the de-signatured and de-multipled data is decomposed into the 
primary PP, PS, SP, and SS wave responses that would be recorded if pure P- and S- 
wave sources and pure P- and S- wave receivers had been used in the seismic survey. 

In one embodiment, Step S6 has two principal operations. Initially, the multi- 
component data are Fourier transformed over horizontal receiver co-ordinates. The 
transformed data then undergo receiver side wave fields decomposition, in the wave 
number domain, into up-going pressure and shear wave fields. This is done by 
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operating on the data by a receiver decomposition matrix. The decomposed data then 
undergoes a reverse Fourier transform, to give the up-going pressure and shear wave 
fields measured at the receiver. 

The thus-obtained multi-component data are then Fourier transformed over horizontal 
source co-ordinates. The transformed data then undergoes source side wave field 
decomposition in the wave number domain, into down- going pressure and shear wave 
fields. This is done by operating on the data with a source decomposition matrix. The 
data then undergoes a reverse Fourier transform, to give the down- going pressure and 
shear wave fields emitted at the source. 

The combined elastic wavefield decomposition on the source and receiver sides gives 
seismic data that are equivalent to data that would be obtained in an idealised survey in 
which the overburden response was absent, and in which single component P- or S- 
wave sources and single component P- or S-wave receivers were used. 

A further advantage of the invention is that it removes, or at least significantly 
attenuates, ground-roll noise (in the case of a land-based seismic survey) or Scholte 
waves (in the case of a marine seismic survey). Ground roll noise is often the largest 
source of noise in a land-based seismic survey, so the removal of the effects of ground 
roll noise is highly desirable. 

Step S6 can be omitted if desired. Alternatively, the step S6 of decomposing the data 
into the PP, PS, SP, and SS wavefields may be carried out before the step S5 of 
applying the de-signature and de-multiple operator. If the decomposition into PP, PS, 
SP and SS wavefields is done first, this step would replace the step S2 of decomposing 
the seismic data into up-going and down-going waves is not required. 

If the medium in which the seismic survey is carried out is horizontally layered, the 
Betti designature and elastic demultiple scheme, followed by the elastic source— receiver 
decomposition scheme, is greatly simplified. It may be conveniently implemented in 
the tau-p domain or in the frequency- wave number domain. 
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The method is preferably carried out on common shot gathers, or on common receiver 
gathers when source array variations are negligible. In the latter domain the assumption 
of a layered earth may possibly degrade the results, but in most cases this is not very 
significant. 

Once the seismic data has been corrected for the effects of the overburden and the 
source radiation characteristics using the de-signature and de-multiple method of the 
invention, and optionally has been decomposed into P-wave and/or S-wave constituents, 
it may then be subjected to further processing steps such as, for example, one or more 
conventional processing steps. 

Figure 4 illustrates the effect of the present invention on a land seismic survey in which 
a multi component seismic source 4 and a multi-component seismic receiver 5 are 
disposed on the earth's surface. 

Figure 4(a) illustrates a real seismic survey. When the source 4 is actuated, down-going 
P- and S-waves are emitted. P-waves are denoted by solid lines, and S-waves are 
denoted by broken lines. The arrow in the box denoting the source (labelled "src") 
indicates the direction of source motion generated by the source. The seismic energy 
emitted from the source 4 undergoes reflection within the earth, and the up-going 
reflected waves are incident on the receiver 5. The receiver 5 contains at least two 
orthogonal geophones, and the arrow in the box denoting the receiver (labelled "rec") 
indicates the component of particle motion that is measured by the receiver. When the 
receiver records the wavefield at the earth's surface it records both p-waves and s-waves 
without distinction. 

Figure 4(a) shows four cases — 

(i) the source generates a horizontal particle motion, and the receiver records particle 
motion in that horizontal direction; 
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(ii) the source generates a horizontal particle motion and the receiver records the 
vertical component of the received particle motion; 

(iii) the source generates a vertical particle motion, and the receiver records a 
horizontal component of the received particle motion; and 

(iv) the source generates a vertical particle motion, and the receiver records the vertical 
component of the received particle motion. 

Figure 4(b) illustrates the results of applying the de-signature and elastic de-multiple 
process to the data recorded in Figure 4(a), That is, Figure 4(b) shows the results after 
the seismic data acquired in the seismic surveying arrangement of Figure 4(a) has been 
processed according to steps SI to S5 of Figure 2. The effect of this processing is to 
eliminate the response of the overburden (that is, the response of the earth's surface in 
this case) from the acquired seismic data. The processing also eliminates the effect of 
the radiation characteristics of the seismic source 5. The source is now a point source of 
seismic energy. In principle energy may propagate upwards past the source and 
receiver, since the overburden has been replaced by an infinite homogenous medium, 
and this is indicated in Figure 4(b). 

As noted above, the source 4 emits both P- and S-waves. In Figures 4(a) and 4(b), the 
receiver 5 records both S-waves and P -waves without distinction. 

Figure 4(c) illustrates the results of applying wave fields decomposition at the receiver 
side to the seismic data of the idealised hypothetical seismic surveying arrangement 
shown in Figure 4(b). This has the effect of replacing the geophones used as the 
receiver in Figures 4(a) and 4(b) by a simulated pure P-wave detector (Figures 4(c)(i) 
and 4(c)(iii)) 5 or by a simulated pure S-wave detector (Figures 4(c)(ii) and 4(c)(iv)). 

Finally, Figure 4(d) shows the results of further applying a source side wavefield 
decomposition. This has the effect of replacing the seismic source used in Figures 4(a) 
to 4(d), which emits both P- waves and S-waves, by a simulated pure P-wave source 
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(Figures 4(c)(i) and 4(c)(ii))or by a simulated pure S-wave source (Figures 4(c)(ii) and 

4(c)(iv)). 

Figure 4(d) illustrates the results of applying steps SI to S 6 of the embodiment shown in 
Figure 2. As can be seen, the method has eliminated the overburden response, 
eliminated the radiation signature of the source, and has decomposed the acquired data 
into PP, PS, SP, and SS events- 
Figure 5 illustrates the results of applying a method of the present invention to an ocean 
bottom seismic survey. In the surveying arrangement shown in Figure 5(a) a seismic 
source 4 and a seismic receiver 5 are each disposed on the sea-bed. The arrow in the 
box denoting the source indicates the direction of source motion generated by the 
source. The receiver 5 contains at least two orthogonal geophones, and the arrow in the 
box denoting the receiver indicates the component of particle motion that is measured 
by the geophone. (A receiver used in a practical ocean bottom seismic survey may also 
contain a hydrophone, but this is not shown in Figure 5.) The four combinations' of 
horizontal and vertical direction of source motion and horizontal and vertical 
component of particle motion measured at the geophone shown in Figures 5(a)(i) to 
5(a)(iv) correspond to the four cases shown in Figures 4(a)(i) to 4(a)(iv). 

As can be seen in Figure 5(a), the source emits down-going S-waves and P-waves, 
which undergo reflection within the earth so as to be incident upon the receiver 5. 
Figure 5(a) also illustrates one possible water layer multiple path for seismic energy, in 
which upgoing P-waves are reflected at the water surface so as to be incident on the 
receiver 5. 

The geophones in the receiver 5 measures the incident S-waves and P-waves without 
distinction. 

Figure 5(b) illustrates the effect of applying a de-signature and elastic de-multiple 
method of the present invention to the seismic data acquired in the seismic surveying 
arrangement shown in Figure 5(a). As can be seen, the water layer multiple paths are 
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eliminated, since the invention eliminates the response of the overburden (that is, the 
water layer in this case) from the acquired seismic data. The method of the invention 
also eliminates the effect of the radiation signature of the source, and the source in the 
idealised seismic surveying arrangement of Figure 5(a) is a point source of force. 

Figure 5(b) illustrates the effects of applying steps SI to S5 shown in Figure 1 to the 
seismic data acquired in the seismic surveying arrangement of Figure 5(a). The effect 
of applying further steps of the invention - that is, applying the receiver side 
decomposition and the source side decomposition would provide identical results to 
those shown in Figures 4(c)(i) to 4(d)(iv). 

Figure 6(a) and 6(b) illustrate the results of applying the present invention to a land 
seismic survey which the receiver and source are deployed within the earth. In the 
seismic surveying arrangement shown in Figure 6(a) the seismic source 4 and the 
seismic receiver 5 are both deployed within the earth, with the receiver being at a 
greater depth than the source. Figure 6(a) illustrates the primary path of seismic energy 
from the source to the receiver, in which down-going S-waves and P-waves axe 
reflected by a target reflector 3 so as to be incident on the receiver 5. Figure 6(a) also 
shows unwanted multiple reflection paths in which up-going S-waves and P-waves 
undergo reflection in layers disposed above the source so as to be incident on the 
receiver. 

In Figure 6(a) the arrows in the box denoting the source 4 indicate the direction of 
source motion generated by the source. The arrows in the box indicating the receiver 
indicate the component of particle motion that is measured by a geophone contained in 
the receiver. As in Figures 4 and 5, the four combinations of horizontal and vertical 
source motion and horizontal and vertical geophone orientation are shown in Figure 
6(a)(i) to 6(a)(iv). 

Figure 6(b) illustrates the effect of applying the present invention to the seismic data 
acquired in the seismic survey arrangement of Figure 6(a). As in Figures 4 and 5, the 
effect of applying a de-signature and de-multiple method of the invention to the data 
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obtained in the seismic survey arrangement of Figure 6(a) is to eliminate the effect of 
the overburden from the acquired seismic data. The effect of the radiation signature of 
the seismic source 5 is also eliminated. Thus, Figure 6(b) illustrates an idealised 
seismic survey in which the overburden response is eliminated (corresponding to an 
homogenous overburden), and in which the seismic source 4 is a point source of force. 
Figure 6(b) illustrates the results of applying steps SI to S5 of Figure 2 to the seismic 
data acquired in the real seismic surveying arrangement of Figure 6(a). 

The effect of applying a receiver side decomposition and a source decomposition would 
correspond to the results shown in Figures 4(c)(i) to 4(d)(iv) respectively. 

Figures 7(a) to 1 1 (d) illustrate results of applying the present invention to simulated 
seismic data. 

Figures 7(a) to 7(d) shows raw simulated seismic data. This data was simulated for a 
land seismic survey in which the source and receiver are both located on the earth's 
surface. The seismic data was simulated using a reflectivity code of the type proposed 
by B. L. Kennett in "Seismic Wave Propagation in Stratefied Media" Cambridge 
University Press (1993). The model used to simulate the seismic data consists of a flat 
free surface on which the source and receiver are located. Below the free surface is a 
300m thick layer, having a flat lower interface. This layer is an isotropic elastic layer, 
with a P-wave velocity of 1800m/s 5 and s-wave velocity of 800m/s, and a density of 
200kg/m . Below this layer is an isotropic elastic infinite half space, having a p-wave 
velocity of 2500m/s > an S-wave velocity of 1500m/s, and a density of 2200kg/m 3 . Since 
the model is horizontally layered, it is necessary to consider a vertical source motion 
and one horizontal source motion, and to consider a geophone that records either a 
vertical component or a horizontal component of the particle motion. Thus, it is only 
necessary to consider a 2C x 2C (4C) case, rather than the full 3C x 3C (9C) case. 

In Figures 7(a) to 10(d), the notation V z and V x indicate vertical particle motion and 
horizontal particle motion at the receiver, respectfully. The terms F z and F x indicate that 
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the geophone records the vertical component and a horizontal component of the particle 
motion, respectively. 

Figure 8(a) to 8(d) illustrates the up-going constituent of the raw data of Figures 7(a) to 
7(d) respectively. Figures 9(a) to 9(d) show the down- going constituent of the raw 
seismic data of Figures 7(a) to 7(d) respectfully. Figures 8(a) to 9(d) were obtained by 
applying equations (A-28), (A-29), (A-32) and (A-33) below. The effect of the free 
surface is still present in Figures 8(a) to 9(d). 

Figures 10(a) to 10(d) illustrate the effect on the seismic data of Figures 7(a) to 7(d) of 
removing the effects of the free surface and the source signature. This is done by 
applying a de-signature and de-multiple method of the present invention. That is, 
Figures 10(a) to 10(d) shows the effect of applying steps SI to S5 of Figure 2 to the raw 
seismic data shown in Figures 7(a) to 7(d), respectively. It can be seen that seismic 
events in Figures 7(a) to 7(d) that relate to multiple reflections have been eliminated 
from the seismic data shown in Figures 10(a) to 10(d). 

The events in the seismic data of Figures 10(a) to 10(d) include all PP, PS, SP and SS 
events, as well as related head-waves. 

Figures 1 1 (a) to 11 (d) illustrate the effect of applying a receiver-side decomposition and 
a source-side decomposition to the data shown in Figures 10(a) to 10(d), respectively. 
The effect of this is to decompose the seismic data into PP 3 SP, PS and SS events. The 
seismic data shown in each of Figures 11(a) to 1 1(d) includes only one event, together 
with the related head-waves. 

The present invention has been described above with reference to a combined de- 
multiple & de-signature method. In an alternative embodiment, the invention provides 
a de-signature method only. In this alternative embodiment, the operator is computed 
from the down-going constituent of the initial arrival of seismic energy at the receiver. 
The initial arrival constitutes seismic energy that has travelled direct from the source to 
the receiver without undergoing any reflection, for example along the path 6 shown in 
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Figure 1(a). The incident part of the wavefield may be determined by, for example, the 
method proposed for acoustic or elastic earth by A. Weglein and E.G. Secrest in 
"Wavelet estimation for a multi-dimensional acoustic or elastic earth", Geophysics Vol. 
55,pp902-913 (1990). 

Furthermore, if the signature/wavelet of the multicomponent source is known or can be 
estimated, the present invention may also be used to provide a de-multiple method only. 
The principal steps of such a method are to (a) obtain a combined de-signature and de- 
multiple operator according to a method of the invention, (b) obtain a de-signature 
operator from the known or estimated radiation characteristics of the seismic source and 
(c) obtain a de-multiple operator that will combine with the de-signature operator 
obtained at (b) to produce the combined de-signature and de-multiple operator obtained 
at (a). 

A detailed mathematical description of one embodiment of the present invention will 
now be given. 



Consider a volume V enclosed by the surface S = 2) + S R with outward-pointing 
normal vector n , as depicted in Figures 3(a) to 3(c). £ is a horizontal plane surface 
located at depth z~ infinitesimally above the multicomponent receiver x r 14 located at 

depth level z r . To simplify the analysis, it is assumed that the solid medium is 

homogeneous and isotropic at depth z r and in a infinitesimally thin region below depth 
2^ The invention does not require knowledge of the properties of the overburden (at 
depths z <z r ) and the subsurface (at depths z > z r ), and the overburden and sub- 
surface may both be arbitrarily inhomogeneous, anisotropic and anelastic media. When 
seismic data is acquired close to a physical surface, the receiver(s) is/are always 
infinitesimally below the physical surface. Physical surfaces typically are the surface of 
the earth (which is assumed to be a solid/vacuum surface), the ocean bottom (which is 
assumed to be a solid/fluid boundary), or another interface in the earth (assumed to be a 
solid/solid surface). S R is a hemispherical surface of radius R . The Cartesian co- 
ordinates will be denoted by x = (<f, x 2 ) , where £ = (x l , x 2 ) . For notational 
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convenience, x 3 = z will also be used. The x 3 -axis, which is positive downwards, is 
parallel to n . The x 1? x 2 -axes are in the 2 plane. 

Initially, an integral relationship is established between the multicomponent source, 
multicomponent receiver data acquired in a real seismic survey - these data contain the 
scattering response of the medium above the receivers - and the desired multicomponent 
source, multicomponent receiver data that would be obtained in an idealised seismic 
survey with that scattering response absent. 

The seismic source used in the survey is assumed to separately generate three 
orthogonal source motions. The desired multicomponent data are those data that would 
be recorded in an idealised, hypothetical multicomponent seismic survey using three 
orthogonal point forces acting separately, with a desired, equal signature when the 
medium above the receivers is homogeneous and extends upwards to infinity, and has 
parameters equal to those of the medium at the receiver depth level. The overburden in 
the idealised survey thus is an elastic, isotropic halfspace. The geology below the 
receiver level is naturally the same in the physical and hypothetical seismic surveys. 





State P 


State H 


State H 


Force component 


AM 


aS(x - x r )S in 


aS[x - x r y> i} 


Traction source 


0 


0 


0 


Particle velocity 


V„(* r |* f ) 


V mn X r ) 




Vertical traction 


S mn {Xr\ X s ) 


S mn X r) 





Table 1: Notation for sources and fields entering the elastodynamic equations 



related to states P, H and H depicted in Figures 3(a) to 3(c), respectively. 

In the physical survey, the seismic source is unidirectional, either a vibrator or a force. 
A vibrator source consists of one vertical vibrator (source of tensile stress) and two 
horizontally oriented vibrators (sources of shearing stress). The vibrator, located at 
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position x s needs not be an ideal point source. Its source signature need not be known. 

In the hypothetical surveys H and H , the sources are point sources of force in the n 
and m direction, respectively, with source signature a . 

Consider the seismic survey with configuration as depicted in Figure 3(a), which will be 
referred to as a "physical" seismic survey for convenience. The recorded rn:th 
component of particle velocity vector at receiver location x r , infinitesimally below 2 , 

due to a source acting in direction n at centre coordinate x s with unknown source 
strength and radiation pattern is denoted by v mn . Likewise, the m:th component of 
vertical traction vector is denoted by s mn . The source and field variables for the physical 

seismic survey, which is denoted for short "state P are listed in Table 1 . Note that 
the surface 2 may or may not coincide with a physical surface. In the case that 
recording takes place just below the surface of the earth, or just below the ocean bottom, 
the surface 2 obviously coincides with a physical surface. 

The desired wavefields, v mn and ? mfl , that it is proposed to solve for are the responses of 

the medium from three orthogonal point forces with desired signature a when the 
medium above the receiver level is a halfspace as shown in Figure 3(b). In this case, the 
surface 2 is a non-physical boundary. The desired particle velocity and vertical 
traction vector responses are recorded at location x r infinitesimally below 2 for the 

point forces located at x' r on 2 . The source and field variables for this hypothetical 
seismic experiment which is denoted, for short, "state H", are listed in Table 1. 

To establish the integral relationship between the physical state P and hypothetical 

state H seismic experiments the hypothetical "state H" is introduced. This is a further 
idealised seismic survey, and corresponds to the idealised seismic survey if if Figure 
3(b) except that the source 5 and the receiver 6 are interchanged with one another. 

Thus, the wavefields v mn and s mr , in the hypothetical "state H" are the reciprocal 

wavefields to the ones in state H , and obey the reciprocity relation 
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v nm (K\ x r)=v mn (x r \x- r ), 0) 

&nm( X ',\ X r)=Sw(xr\x~,)> ( 2 ) 



Thus, v w>i and S mn are responses at location x' r on surface E due to a point force, with 

signature a , acting in direction at location je r infinitesimally below £ as depicted 

in Figure 3(c). The surface E is, as in the desired idealised state H , an artificial, non- 
physical boundary. Below, we shall use the fact that seismic data acquired in a survey 

in the hypothetical state H consist of upgoing events only scattered from the subsurface 
below E . In addition, the direct wavemodes from the sources to the receivers are 
upgoing events since the sources are infinitesimally below the receivers. The source 

and field variables for state H are listed in Table 1 . 



Next Betti's reciprocity theorem is applied in the volume V enclosed by the surface 
E + S R ^ , where is located at infinity. Similarly to, for instance L. Amundsen et 

al. in "Elimination of free-surface related multiples without the need of the source 
wavelet", Geophysics, Vol 66, pages 327-341 (2000), a frequency-domain integral 

equation describing the relationship between state P and state H is obtained: 
av WM (* r |xJ= . (3) 



The surface S R does not contribute to the integral as R — > oo (this is the radiation 

condition, see Y.H. Pao and V. Varatharajulu in "Huygen's Principle, radiation 
conditions and integral formulae for the scattering of elastic waves", J, Acoust. Soc. Am. 
Vol. 59, ppl361-1371 (1976)). Equation (3) can be simplified by identifying proper 
boundary conditions for E . In the physical state P , v in and s^ are sums of up-going 

and down-going waves: 



(4) 
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in in in 



(5) 



whereas in the hypothetical state H , v. and f . are purely up-going fields: 



v. = v (tt) ; v (d) =0 

in in * in 



s. = s w ; = 

in in ? m 



(6) 
(V) 



These boundary conditions are most conveniently introduced into equation (3) by 
analysing the problem in the horizontal wavenumber domain, where up -going and 
down-going waves and their relation to vertical traction and particle velocity vectors are 
explicitly known. Making use of Parceval's identity : 



OO -t 00 



\dKF{- k)h(k) , (8) 



-co v y -co 



where wavenumber vector k = {k x ,k 2 } is conjugate to % = (x l ,x 2 ), equation (3) yields 



a 



v mn (x, \x s ) = — l — \dK:[s im (k, z~ r \x r Y in (- k z r \x s )- V im («; z' r \x r )s in (- k; z r |x, )] . 



(9) 



Going from summation convention to vector notation, equation (9) is written 



av mn (x, \x s )= -^—^ \dK[S T m [k, z~ r \x r y» (- k, z, \x s )- V T m {fc 3 z~ \x r )s n (- k z, \x s )] 3 (10) 



where S T = (S l , S 2 , 5 3 ) and V T - (v x , V 2 , V l ) are wavenumber domain vertical traction 
vector and particle velocity vector, respectively, and the superscript T means transpose. 

As shown in the Appendix, since the hypothetical state H fields S and V consist of 
up-going wave modes only, they are related as 
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s(k) = L sv (k)P(k:) , (11) 



where L sv is a 3x3 matrix depending on the local medium parameters p , a , jB along 

the receiver spread. The elements of L sv are given in equation (A- 19) in the appendix. 
Inserting equation (11) into equation (10) gives 



QV ™ (*, \ X s )= \ dK K (^ Z 'r\ X r ^ (~ K Z, \x, ) , (12) 



where 



R m =L sy V-S = G 1 V (D) 



(13) 



is interpreted as the overburden response due to scattering in the medium above the 
receiver level. Furthermore, the vector 

V m =[vt D \Vt D \Vj D)1 [ (14) 



contains the elements of the down-going wavemodes on each of the particle velocity 
components V x , V 2 , and F 3 . Generally, for every shot location, V m is calculated in 

the slowness or wavenumber domain from the particle velocity vector and the vertical 
traction vector according to equations (A-29), (A-31) and (A-33). The elements of the 
resulting down-going reflection response are: 



*r - J 2 ? k -pifi 2 Q a -iJrr+P^ 2 ((! a -4 fi >r} as) 

p +9 a 9 fi 



Ri D) = 2 2p W? 2 - 9, >r + k - pIP t - 9 J ^ (D) } (16) 
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^3 ~ 2 . K 3 



(17) 



Re-introducing the summation convention, we obtain 



a v mn (x r \x s )= \dK V im (a; z; \x r (- k 3 z r \x s ) . (18) 



Making use of Paxceval's identity yields in the space domain 

2-v M „(*,|*J= [dS(z)v lm {x\x r yw{x\x,), (19) 
where is the Fourier transform of . 

Equation (19) gives an integral relationship between the field v nm in the hypothetical 

state H and the recorded field v rm in the state P . Note that no information, except 

location, about the physical source and its radiation characteristics, and no information 
of the properties of the physical overburden above the surface S or the physical 
subsurface below the receiver level have been used to derive integral equations 
(19) for v nm . The hypothetical state H fields v nm are directly determined from 
reciprocity. 

Observe that the inverse of rf n d) in equation (19) can be interpreted as a 

multidimensional operator that acts as (i) a deterministic designature operator which 
removes the physical source characteristics, and (ii) a deterministic multiple attenuation 
operator that eliminates the overburden response from the physical data. 

Eliminating the incident wavefield in the hypothetical state 
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The desired field v nm in the hypothetical experiment can be split into an incident wave 

field v^ c) propagating upwards from the source to the receiver, and the wavefield v£f 
scattered upwards from the subsurface: 



nm nm nm 



where G is the Green's tensor, equivalent to the tensor derived in equation (A-39) 
below. Observing that 

(v™ J (k, z- \x r )R w (- k, z r \x, )= SG T (k^ 1 (- K y w (- «r> X p(- i K £ ) 

= aV m (- k )exp(- iK ) 



and 



we obtain: 



(20) 



The incident wave field, which propagates in a homogeneous medium and thus is 
independent of source location, can be shown to be given by: 

y (inc) k k ) = a G exp(- i K ■ £ ) , (21) 



(22) 



v£ } (*, K ) = -^rr J** «p(- i* • £ {- k, z, \x s ) , (23) 



dlv^ r \x t )=-^jdK^{^z;\x r ^(-^z r \x t ) . (24) 



Making use of Parceval's identity yields in the space domain: 

«v£> (x, \x t )=[ dS{$) (x|x r K D) (*|*. ) • (25) 
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Equation (25) gives the sought-after integral relationship between the scattered field 
v^ c) in the hypothetical state H survey and the up-going and down-going fields 

and in the state P . The scattered fields v^ c) in the hypothetical state H are 
determined from the reciprocity relation. 

Equation (25) is a Fredholm integral equation of the first kind for the desired scattered 
fields, leading to a system of equations that can be solved for v^ c) by keeping the 
receiver co-ordinate fixed while varying the source co-ordinate. 



Wavenumber domain solution 



Fourier transforming equation (25) over source coordinates £ s and receiver coordinates 
% r yields 



dc r , z r \k s , z s )= ^JL_ [dK V™ fa z~ r \K r , z r > ™ (- k, z, \k s , z, ) . (26) 



Equation (26) leads to a system of equations that can be solved for by keeping the 

wavenumber conjugate to the receiver co-ordinate fixed while varying the wavenumber 
conjugate to the source co-ordinate. The coupling between positive wavenumbers in the 
down-going overburden response field with negative wavenumbers in the desired field 
(and vice versa) reflects the autocorrelation process between the two fields. 

Decomposition into primary PP, PS, SP, and SS wave responses 

Following the Betti designature and elastic demultiple step is an elastic wavefield 
decomposition step that decomposes the multicomponent source, multicomponent 

receiver Betti designatured and elastic demultipled data v (iC) (jc*|jc* ) given in equation 
(25) into primary PP, PS ? SP, and SS wave responses that would be recorded from 
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pressure wave and shear wave sources and receivers. The receiver and source co- 
ordinates for the hypothetical data are denoted by x h r and x* , respectively, where x h r is 

a point on the surface S and x* = x r is infinitesimally below. Let matrix 0 contain 
up-going potentials at the receiver from sources of force in direction i , 





[ u„ 


o P2 


u P3 ) 






U SV1 


U SV2 


U SV3 


(27) 






^ SH2 


U SH3 y 





Let U be the matrix containing up-going P, SV, and SH wave potentials at the receiver 
excited from P, SV, and SH wave sources, 





r U PP 


U psv 


u PSH 


u = 


U SVP 


U svsv 


U SVSH 




^SHP 


U SHSV 


U SHSH J 



(28) 



This elastic source-receiver wavefield decomposition may now be divided into two 
computational operations. First, the multicomponent data v ( "° are Fourier transformed 

over horizontal receiver co-ordinates £ * . In the wavenumber domain K r , receiver side 
wavefield decomposition into u-pgoing pressure and shear wavefields is achieved by 
left-multiplying V (sc) by a receiver decomposition matrix R , 

0(K„z h \x h s )=R(K r )V <*> (k, ,z h r \x h s ), (29) 



where 



WO 03/036331 



PCT/IB02/05056 



M&C Folio No P5 1685 WO 



35 



f 



R 



1 



p +q*<ip 



V 



Pi 

PlQa 



Pi 

P 2 <la 



P P 
P P 



~P 

0 



J 



(30) 



with p~K r I co . An inverse Fourier transform gives upgoing pressure and shear 
wavefields measured at the receiver u[x h r \x h s ). Second, the multicomponent data 

u[x h r ) are Fourier transformed over horizontal source co-ordinates £ * . In the 

wavenumber domain k s , source side wavefield decomposition into down-going 

pressure and shear source wavefields is achieved by right-multiplying the data by a 
source decomposition operator 5 , 



u(x h \K s> Z h s )=u(x h \K s ,Z>;)s{K s ) , 



(31) 



where 



f 



S = 



-ia> l p 

p 2 +q a <i0 





PiVa 






P 




Plla 


Pi 






P 




~P 



PlQa 



P 
PlQa 



(32) 



0 



V 



with p — k s I co . An inverse Fourier transform gives data U equivalent to data from a 

hypothetical survey with overburden absent, with single component pressure and shear 
wave sources, and single component pressure and shear wave receivers. 

In {jc r , jc s ) space, the source-receiver decomposition is 
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x: r ,z r k s ,z s 



(33) 



Observe that the decomposition operators are simply matrices which are multiplied by 
the data. The required parameters for the decomposition are the density, and P-wave and 
S-wave velocities at the receiver depth level. 

Betti decon volution and source-receiver wavefield decomposition: Horizontally 
layered medium ("1.5D medium") 

In a horizontally layered medium, the response is laterally shift-invariant with respect to 
horizontal source location £ 5 , and the expressions therefore simplify significantly. For 

this reason, a horizontally layered medium 1 s sometimes referred to as a "1 .5D 
medium". Assuming a horizontally layered medium often yields satisfactory results for 
real data and realistic Earth structures. 

Consider equation (24), where we may set g s = 0 and 



By Fourier transforming equation (24) with respect to £ r and interchanging integrals, 
we find 



The last integral is recognised as the Dirac delta function. Performing the integration 
over wavenumbers, using the Dirac delta function property, and renaming K r by k , we 

obtain 




(34) 




(35) 
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2VZ,(k, Z ,\Z s =0,z s )=V£i-K,z± = 0,z r )^(K,z;\Z i = 0,z s ) 



(36) 



Writing this result in terms of matrices, yields 



( vOc) 
y 11 



(sc) 
12 



T?(™) 



yi.sc) 

r 22 
K 23 



yv 

r 31 
r 32 
*33 



\ 



= a 



V 



yW) y 

V 21 V ' 



(U) 
22 

yw 

v 32 



yW\ 

r 13 



r 23 

(£/) 



33 



i?^ R 



iV 21 



2? 



31 



(Z0 

12 
1V 22 
iY 32 



R 



3 

CD) 
23 

33 y 



(37) 



Using the field properties 




it is evident that the scattered part of the desired field is obtained by generalised spectral 
deconvolution: 



Y ii 

v 12 



M3 K : 



' v 21 

rr{sc) 
V 22 

(sc) 
23 



r 31 



32 



*33 



yw 

Y 21 



*12 

22 
32 



r n 

V 



13 

(CO 

23 



r 33 y 



V, 



\ 



^11 



JV 21 

K 31 



12 
iV 22 
lX 32 



R 

R 
R 



13 

CD) 
23 



-i 



33 



(■.*r|*r) 

(39) 



It follows from equation (39) that the components of the desired field are obtained by 
deterministic spectral deconvolution between the field itself and the reflection response 
of the overburden containing the down-going part of the particle velocity vector. 

The next step is to decompose the Betti deconvolved data into primary PP, PS, SP 5 and 
SS wave responses. This is carried out as described in the section "Decomposition into 
primary PP, PS, SP ? and SS wave responses". 

Observe that the 1.5D Betti deconvolution scheme followed by source-receiver 
decomposition may be implemented as r - p or frequency-wavenumber domain 
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algorithms. In the t — p domain, a joint designature, multiple attenuation and source- 
receiver decomposition process is performed for each p -trace. In the frequency- 

wavenumber domain, the process is performed for each combination of frequency and 
wavenumber. 

Applications for the present invention 

As noted above, the present invention can be applied to a number of different seismic 
survey configurations. Three survey scenarios of particular interest will be discussed 
briefly. 

1. Land surface seismic acquisition 

The invention maybe applied to a 3Cx3C (or 9C) land seismic survey where three 
orthogonal source motions (two horizontal and one vertical) are generated separately, 
and where the seismic wavefield is recorded by three orthogonal geophones, measuring 
the two horizontal and the vertical components of the particle velocity vector. 

For measurements located at the free surface, the up/down separation simplifies 
considerably: all tractions S l , S 2 , and S 3 vanish so that the traction terms vanish in 

equations (A-28) to (A-33) below. 

2. Seabed seismic acquisition 

The invention may be applied to seismic data acquired in a 3Cx4C (or 12C) seabed 
seismic experiment where three orthogonal source motions (two horizontal and one 
vertical) are generated separately, either on the ocean bottom or in the water column, 
and where the seismic wavefield is recorded at a receiver location by three orthogonal 
geophones deployed just below the sea floor and one hydrophone deployed just above 
the sea floor. 
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For measurements located at the fluid/solid interface of the seabed, the up/down 
separation simplifies considerably: tractions iS' 1 and S 2 vanish so that the corresponding 

traction terms vanish in equations (A-28) to (A-33) below. The hydrophone 
measurement on the seafloor is the negative of S 3 . 

It is worthwhile to note that sources that are equivalent to sources of force can be 
generated in the water by using, for instance, a conventional airgun source. For 
instance, the response corresponding to a vertical point-force can be obtained by 
acquiring two records with the same source at the same location but at slightly different 
depths. A finite-difference approximation then allows us to compute the response due 
to a pressure gradient source in the vertical direction. The equation of motion 
(Newton's second law) shows that this is equivalent to a vertical point force. The same 
is true for the two horizontal directions, as well. 

3. Borehole seismic acquisition 

The 3C x 6C (or 18C) borehole seismic experiment where the three orthogonal source 
motions are generated separately, either on land, on the sea floor, or in the water 
column, for which the three components of particle velocity and the three components 
of vertical traction are known. 

In this surveying arrangement, there are no boundary conditions that simplify the 
up/down separation step. However, by adding hydrophones in the receiver package, 
additional constraints can be obtained (a hydrophone measures divergence of particle 
motion). 

Wavenumber domain basic relationships between the particle velocity and vertical 
traction vectors and upgoing and downgoing wave vectors 

This section considers a horizontally layered elastic earth. In a source-free region 
elastic wave propagation is described by the equation of motion and Hooke's law (the 
elastic constitutive relation). These can be written as a system of first-order ordinary 
differential equations of the form 
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d,B = -io)AB 



where the field vector B is defined 



B = {V T ,S T J 

with particle velocity vector V T = [V, ,V 2 ,V J \ and vertical traction vector 
S 7 — [l?! , S 2 , S3 J . The elastic system matrix A has the form 



(A-l) 



(A-2) 



r 



A = 



0 



0 



0 



0 



A 



A 



A + 2ju 

-0p l p 2 
0 



A + 2/U 



p-9p 2 -jjp a p g 
0 
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Pi 






M 
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P2 






M 
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0 
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0 


0 


0 


0 


0 


p 


Pi 


Pi 



0 
0 

1 



A + 2/u 
A 

P, 



A + 2ju 1 
A 

p 



A + 2{i 
0 



where 6 = p(3A+2 //)/( A +2 /J) and p 2 = p\ + pi . 



(A-3) 



For notational convenience, the explicit dependence of different quantities on 
frequency, wavenumber, depth, etc., is omitted. For instance, the particle velocity vector 

recorded at depth x 3 , V^x^eyjXy) > due to a point source at location X s is in the . 

wavenumber domain denoted V or V{x 3 ) with the understanding 

When required we will show dependency on horizontal 
slowness vector p = k/co = (p^p 2 )' 



Up-going and down-going waves 



For the decomposition of the elastic field into up- and down-going waves in a layered 
earth it is necessary to find the eigenvalues and eigenvectors of the system matrix A for 
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given wavenumbers and frequencies. The field vector B can be decomposed into up- 
(17) and down- going (D) waves 

JV = [U T ,D T ] T (A-4) 



where U T = [f/ p> U sv , U SH ] and D T = [D P , D sv , D SH ], by the linear transformation 
B = LW 

where L is the local eigenvector matrix of A (i.e., each column of L is an 
eigenvector). Equation (A-5) describes composition of the wavefield B from its 
up going and downgoing constituents. 

Given the inverse eigenvector matrix L" 1 y the up- and downgoing waves can be 
computed by evaluating 



(A-5) 



W = L 1 B . 



(A-6) 



Equation (A-6) describes decomposition of the wavefield B into upgoing and 
downgoing P- and S-waves. After some straightforward, but tedious calculations, the 
composition matrix is obtained 

L = L(p)=( Lyu ^ ~ L w(~Pf\ 
\ L su(p) L su(~P)j ' 

and the decomposition matrix 

f lT su{p) lT vu(p) 



(A-7) 



where 
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1 



P2 Qfi 
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Pi 1 



4ml p\p 



a 
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V 



[m~fi 



p [p<[fi 
o 



1 



72 



-2fjp x 
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In a source- free homogeneous solid upgoing and downgoing waves satisfy the 
differential equations: 

d 3 U P =-i<aq a U p , 

d i u sH =-i03q p U SH , 
djDj, = icoq a D p , 
d } D sr =ia>q p D sr , 

d 3 D SH = icoq p D SH . 

Traction-particle velocity vector relationship when down-going waves vanish 
Using the relation (A-5) with 



(A-9) 



(A- 10) 



(A-ll) 



(A- 12) 



(A- 13) 



(A- 14) 



(A- 15) 



(A- 16) 



Z> = 0 



one obtains a simple relationship between the vertical traction vector and the particle 
velocity vector, 



WO 03/036331 



PCT/IB02/05056 



M&C Folio NoP51685WO 



43 



S = L SV V, 



where 



J\ A. 



T — T 7~ l 



(A- 17) 



(A- 18) 



r 



p +q«<ifi 



v 



9a -plP 2 (q a -9b) 
PiPiP 2 {<i a -qp) 

p l l-2/3 2 (p 2 +q a q /3 J 



PlPz0 2 (<la Pi 
9a ~P 2 lP%a ~9fi) Pi 

p 2 l-2p 2 (p 2 +q a qpl 



l-2j3 2 (p 2 + q a q p 
1 - 2 j8 2 {p 2 +q a q /3 

9p 



(A- 19) 



Upgoing and downgoing waves scaled as particle velocity and vertical traction 

Upgoing and downgoing P, SV, and SH waves are not uniquely defined. We may scale 
the components so that they have dimensions of particle velocity or traction in such a 
way that their sum gives a component of particle velocity or a component of vertical 

traction. Let V^ Up) denote the upgoing P-waves on V. , VS Usv) denote the upgoing SV- 
waves on V { , etc. Then, 



y — y(Pp) y( u sv) _|_ y( u sx) _j_ yi D p) _j_ y( D sr) y( D sft) 



i i 



with similar equation for vertical traction. Furthermore, it is possible to sum the upgoing 
components into a total upgoing component, and likewise for the downgoing 

components. V^ U) is defined as the sum of upgoing waves on V i , 



y(U) _ y{U P ) + y{U SY ) + y{U SH ) ^ 



and Vf D) as the sum of downgoing waves on V i , 



yW = y(Dp) + y&sv) + y&sif) t 



The total upgoing waves on S £ are 
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o(U) _ «{U P ) n(U sy ) r,(Cf w ) 
I i i i ' 

and the total downgoing waves on S t , 



Further, defining the vectors 



y«o=[y<m,VW> t yw>Y (A . 20 ) 

F w =[F r I w ,F I (0) ,F 3 w f (A-21) 
S^=[^,Sf\.?rF (A-22) 
Sf w = [S™ , , > f (A-23) 

it follows that these vectors are related to the originally defined vectors of upgoing and 
downgoing waves, U and D , respectively, as 

V (V) ip)=L yu (p)u{p) (A-24) 
V « (/>) = (- p)2)(p) (A-25) 

S ra (/>) = £ W (/>M/>) (A-26) 
5 w (p) = X JP (-/,)D(p) (A-27) 

K w and enter the integral relationship between the physical experiment and the 
hypothetical experiment as known field vectors. 



Explicitly, the total upgoing and downgoing waves on V l are given from V x , V 3 , and 
as 

^ -£-l-2/3% 2 + T i -(p a + (A-28) 
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K D) = ^ + ^ |I - 2^ (p 2 -h q a q p )> 3 " (p 2 + q a q, X ■ ( A " 29 > 

Note that total up-going and down-going waves on V x do not depend on V 2 , 5 2 , and 
1S3 . Further, total up-going and down-going waves on V 2 are given from V 2 , V 3 , and 
as 

VP =^~{ 1 -l-2jB 2 (p 2 +qaq0 f 3+ -^—( P 2 +qaq ,y i , (A-30) 

V 2 (D) =\v 2 +i^-[-2/3*(p> + q a q^--~ 1 -(p 1 +q a q^> (A-31) 

and are thus independent of , iS 2 , and 5 3 . Total up-going and down- going waves on 
V 3 are given from V x , V 2 , V % , and S 3 as 

^ {U) = ± + — — 2 + <7 a 9 , X + f 1 - 2 >S 2 (P 2 + ^ 1a ^ + P 2 V2 1 (A- 32) 

2 2pq fi 2q 0 

Vr =\v, --^—(p 2 + --Ml-2/? 2 ^ 2 +^]Cp 1 K 1 + J p a K 2 ). (A-33) 

2 2p^ 2^ 

Derivation of matrix relations 

Equation (A- 18) and symmetry relations of (A- 19) yields 

[i S¥ v -s\p)= l£j ( p )l t su {- P )v(p) - s(p) 

= L-^(-piL T su (-py(j,)-L T ru (-p)s(p)] ' 



(A-34) 



The last term in the square brackets is identified by use of equations (A-6) and (A-8) as 
the downgoing wave vector calculated from V and S according to 

D{p>L T su yy{p)-I^j{-p)S{jp) . (A-35) 



In the previous subsection, it was shown that D was related to the down-going particle 
velocity wave vector V m by equations (A-28) and (A-29). Inverting this equation gives 
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D{p) = -LA,{- P y m (p). (A-36) 
which is substituted into equation (A-34) to give 

[L sy V-sip) = L^(-p)LUp)V (D) {p) 

= \L m (:P% T ya (-/>)]' V m GO , (A-37) 
= G 1 (py™ (p) 

where 

G 1 =^ vu L T yu )' . (A-38) 
Explicitly, it is found that 



p 2 + q a <ip 



V 0 0 q pj 



(A-39) 



In the theory given above it has been assumed that the sources in the idealised survey 
are point sources of force. However, the invention is not limited to the use of sources in 
the idealised survey which are point sources of force, although the detailed theory given 
above requires slight modification to accommodate sources in the idealised survey that 
are not point sources of force. 

In the detailed description given above the designature/demultiple operator is derived 
using Betti's theorem. The invention is not however limited to deriving the designature/ 
demultiple operator using Betti's theorem. The designature/demultiple operator of the 
invention may be derived from the elastodynamic wave equation or a representation of 
the elastodynamic wave equation. Betti's theorem is an integral representation of the 
elastodynamic wave equation, and examples of other representations of the 
elastodynamic wave equation include the elastic Kirchhoff integral or the elastodynamic 
representation theorem, and the inverse scattering series method. 
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Figure 12 is a schematic block diagram of an apparatus 15 according to the 
present invention. The apparatus is able to carry out a method according to the 
present invention. 

The apparatus 15 comprises a programmable data processor 16 with a program 
memory 17, for instance in the form of a read only memory ROM, storing a 
program for controlling the data processor 17 to perform a method of the 
invention. The system further comprises non- volatile read/write memory 18 for 
storing, for example, any data which must be retained in the absence of power 
supply. A "working" or "scratchpad" memory for the data processor is provided 
by a random access memory (RAM) 19. An input device 20 is provided, for 
instance for receiving user commands and data. An output device 21 is 
provided, for instance for displaying information relating to the progress and 
result of the method. The output device may be, for example, a printer, a visual 
display unit or an output memory. 

Seismic data for processing according to a method of the invention may be 
supplied via the input device 20 or may optionally be provided by a machine- 
readable store 22. 

The program for operating the system and for performing the method described 
hereinbefore is stored in the program memory 17, which may be embodied as a 
semi-conductor memory, for instance of the well-known ROM type. However, 
the program may be stored in any other suitable storage medium, such as 
magnetic data carrier 17a (such as a "floppy disc") or CD-ROM 17b. 
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CLAIMS: 

1 . A method of processing multi-component seismic data, the data having been 
acquired by emitting multi-component seismic energy at a source location; and 
acquiring seismic data at a multi-component seismic receiver located at a greater depth 
than the source location, the method comprising the steps of: 

decomposing the acquired seismic data into up-going constituents and down- 
going constituents; and 

calculating a de-signature and de-multiple operator from the down-going 
constituents of the acquired seismic data and from properties of the medium 
surrounding the receiver. 

2. A method as claimed in claim 1 and comprising the further step of processing 
the acquired seismic data using the de-signature and de-multiple operator thereby to 
attenuate or remove seismic events arising from multiple reflections. 

3. A method as claimed in claim 2 wherein the step of processing the acquired 
seismic data comprises processing the up-going constituents of the acquired seismic 
data using the de-signature and de-multiple operator. 

4 A method as claimed in claim 2 or 3 wherein the step of processing the acquired 
seismic data further comprises selecting a desired seismic signature for the source. 

5. A method as claimed in claim 2, 3 or 4 and comprising the further step of 
decomposing the processed seismic data into P-wave and/or S-wave data. 

6. A method as claimed in claim 1 and comprising the step of decomposing the 
acquired seismic data into P-wave and/or S-wave data, the step of decomposing the 
acquired seismic data into P-wave and/or S-wave data being performed before the step 
of decomposing the seismic data into up-going constituents and down-going 
constituents. 
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7. A method as claimed in claim 5 or 6 wherein the step of decomposing the 
seismic data into P-wave and/or S-wave data comprises decomposing the data at the 
receiver side. 

8. A method as claimed in claim 5, 6 or 7 wherein the step of decomposing the 
seismic data into P-wave and/or S-wave data comprises decomposing the data at the 
source side. 

9. A method of processing multi-component seismic data, the data having been 
acquired by emitting multi-component seismic energy at a source location; and 
acquiring seismic data at a multi-component seismic receiver located at a greater depth 
than the source location, the method comprising the steps of: 

decomposing the acquired seismic data into up-going constituents and down- 
going constituents; and 

calculating a de-signature operator from the initial down-going constituents of 
the acquired seismic data and from properties of the medium surrounding the receiver. 

10. A method of processing multi-component seismic data, the data having been 
acquired by emitting multi-component seismic energy at a source location; and 
acquiring seismic data at a multi-component seismic receiver located at a greater depth 
than the source location, the method comprising the steps of: 

decomposing the acquired seismic data into PP-constituents, PS-constituent, SP- 
constituents and SS-constituents; and 

calculating a de-signature and de-multiple operator from the constituents of the 
acquired seismic data and from properties of the medium surrounding the receiver, 

11. A method as claimed in claim 9 or 10 and comprising the further step of 
processing the acquired seismic data using the de-signature and de-multiple operator 
thereby to attenuate or remove seismic events arising from multiple reflections. 

12. A method of calculating a de-signature and de-multiple operator for multi- 
component seismic data, the data having been acquired by emitting multi-component 
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seismic energy at a source location; and acquiring seismic data at a multi-component 
seismic receiver located at a greater depth than the source location, the method 
comprising the steps of: 

decomposing the acquired seismic data into up-going constituents and down- 
going constituents; and 

calculating a de-signature and de-multiple operator from the down-going 
constituents of the acquired seismic data and from properties of the medium 
surrounding the receiver. 

13. A method of calculating a de-signature operator for multi-component seismic 
data, the data having been acquired by emitting multi-component seismic energy at a 
source location; and acquiring seismic data at a multi-component seismic receiver 
located at a greater depth than the source location, the method comprising the steps of: 

decomposing the acquired seismic data into up-going constituents and down- 
going constituents; and 

calculating a de-signature operator from the initial down-going constituents of 
the acquired seismic data and from properties of the medium surrounding the receiver. 

14. A method of calculating a de-signature and de-multiple operator for multi- 
component seismic data, the data having been acquired by emitting multi-component 
seismic energy at a source location; and acquiring seismic data at a multi-component 
seismic receiver located at a greater depth than the source location, the method 
comprising the steps of: 

decomposing the acquired seismic data into PP-constituents, PS -constituent, SP- 
constituents and SS-constituents; and 

calculating a de-signature and de-multiple operator from the constituents of the 
acquired seismic data and from properties of the medium surrounding the receiver. 

15. A method as claimed in any preceding claim and further comprising the steps of: 
emitting multi-component seismic energy at a source location; and acquiring seismic 
data at a multi-component seismic receiver located at a greater depth than the source 
location. 
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16. An apparatus for processing multi-component seismic data, the data having been 
acquired by emitting multi-component seismic energy at a source location; and 
acquiring seismic data at a multi-component seismic receiver located at a greater depth 
than the source location, the apparatus comprising: 

means for decomposing the acquired seismic data into up-going constituents and 
down- going constituents; and 

means for calculating a de-signature and de-multiple operator from the down- 
going constituents of the acquired seismic data and from properties of the medium 
surrounding the receiver. 

17. An apparatus as claimed in claim 1 6 and further comprising means for 
processing the acquired seismic data using the de-signature and de-multiple operator 
thereby to attenuate or remove seismic events arising from multiple reflections. 

18. An apparatus as claimed in claim 17 and further comprising means for 
processing the up-going constituents of the acquired seismic data using the de-signature 
and de-multiple operator. 

19. An apparatus as claimed in claim 16, 17 or 18 and further comprising means for 
selecting a desired seismic signature for the source. 

20. An apparatus as claimed in claim 17 and further comprising means for 
decomposing the processed seismic data into P-wave and/or S-wave data. 

21. An apparatus as claimed in claim 20 and further comprising means for 
decomposing the seismic data into P-wave and/or S-wave data at the receiver side. 

22. An apparatus as claimed in claim 20 or 21 and further comprising means for 
decomposing the seismic data into P-wave and/or S-wave data at the source side. 

23 . An apparatus for processing multi-component seismic data, the data having been 
acquired by emitting multi-component seismic energy at a source location; and 
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acquiring seismic data at a multi-component seismic receiver located at a greater depth 
than the source location, the apparatus comprising: 

means for decomposing the acquired seismic data into up-going constituents and 
down-going constituents; and 

means for calculating a de-signature operator from the initial down-going 
constituents of the acquired seismic data and from properties of the medium 
surrounding the receiver. 

24. An apparatus for processing multi-component seismic data, the data having been 
acquired by emitting multi-component seismic energy at a source location; and 
acquiring seismic data at a multi-component seismic receiver located at a greater depth 
than the source location, the apparatus comprising: 

means for decomposing the acquired seismic data into PP-constituents, PS- 
constituent, SP-constituents and SS-constituents; and 

means for calculating a de-signature and de-multiple operator from the 
constituents of the acquired seismic data and from properties of the medium 
surrounding the receiver. 

25. A method as claimed in claim 9 or 10 and comprising the further step of 
processing the acquired seismic data using the de-signature and de-multiple operator 
thereby to attenuate or remove seismic events arising from multiple reflections. 

26. An apparatus for calculating a de-signature and de-multiple operator for multi- 
component seismic data, the data having been acquired by emitting multi-component 
seismic energy at a source location; and acquiring seismic data at a multi-component 
seismic receiver located at a greater depth than the source location, the apparatus 
comprising: 

means for decomposing the acquired seismic data into up-going constituents and 
down-going constituents; and 

means for calculating a de-signature and de-multiple operator from the down- 
going constituents of the acquired seismic data and from properties of the medium 
surrounding the receiver. 
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27. An apparatus for calculating a de-signature operator for multi-component 
seismic data, the data having been acquired by emitting multi-component seismic 
energy at a source location; and acquiring seismic data at a multi-component seismic 
receiver located at a greater depth than the source location, the apparatus comprising: 

means for decomposing the acquired seismic data into up-going constituents and 
down- going constituents; and 

means for calculating a de-signature operator from the initial down-going 
constituents of the acquired seismic data and from properties of the medium 
surrounding the receiver. 

28. An apparatus for calculating a de-signature and de-multiple operator for multi- 
component seismic data, the data having been acquired by emitting multi-component 
seismic energy at a source location; and acquiring seismic data at a multi-component 
seismic receiver located at a greater depth than the source location, the apparatus 
comprising: 

means for decomposing the acquired seismic data into PP-constituents, PS- 
constituent, SP-constituents and SS-constituents; and 

means for calculating a de-signature and de-multiple operator from the 
constituents of the acquired seismic data and from properties of the medium 
surrounding the receiver. 

29. An apparatus as claimed in any of claims 16 to 28 and comprising a 
programmable data processor. 

30. A storage medium containing a program for a data processor of an apparatus as 
defined in claim 29. 
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Particle velocities with free-surface present 

Raw data (2C source and 2C receiver) 
generated using a reflectivity code(Kennett,1983). 
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Particle velocities with free-surface absent 
Particle velocities after applying the present invention 

to remove me effects of the free surface. 
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PP,SP,PS and SS wavefields after applying the present invention to decompose 
particle velocities on the receiver and source side into P-wave and 

S-wave receivers and sources. 
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